GLOBAL SOLVABILITY AND BLOW UP FOR THE CONVECTIVE 
CAHN-HILLIARD EQUATIONS WITH CONCAVE POTENTIALS 



A. EDEN, V.K. KALANTAROV AND S.V. ZELIK 

Abstract. We study initial boundary value problems for the convective Cahn-Hilliard 
equation dtu + d*u + ud x u + d x (\u\ p u) = 0. It is well-known that without the convective 
term, the solutions of this equation may blow up in finite time for any p > 0. In contrast 
to that, we show that the presence of the convective term ud x u in the Cahn-Hilliard 
equation prevents blow up at least for < p < |. We also show that the blowing up 
solutions still exist if p is large enough (p > 2). The related equations like Kolmogorov- 
Sivashinsky-Spiegel equation, sixth order convective Cahn-Hilliard equation, are also 
considered. 



1. Introduction 

It is well-known that the solutions of the semilinear heat equations with concave po- 
tentials 

d t u — Au — u\u\ p = 

blow up in finite time if p > and the initial energy is negative, see e.g. [1, 17, 21, 20, 28] 
and references therein. However, it is also established that the presence of the convective 
terms in the semilinear parabolic equation prevents blow up if the nonlinear source term 
is not growing very rapidly. For instance, the solutions of the following convective heat 
equation 

d t u + ud x u — d 2 x u — u\u\ p = (1.1) 

in a bounded interval Q = [-L, L] with Dirichlet boundary conditions exist globally in 
time if p < 1 and the blowing up solutions occur only if p > 1, see [5, 22, 23, 30], see 
also [32] for the results on suppressing the blow up by adding the sufficiently large linear 
convective terms in reaction-diffusion equations. 

The main aim of the present paper is to study the analogous problems for the following 
convective Cahn-Hilliard (CH) equation with concave potentials 

d t u + d 2 x (d 2 x u + u\u\ p ) + ud x u = (1.2) 

in a bounded segment Q = [-L, L] endowed by periodic boundary conditions. 

The long-time behavior of solutions of initial boundary value problems for CH and 
related equations are intensively studied by many authors, see [2, 4, 10, 16, 14, 26, 31] 
and references therein. For instance, the existence of the blowing up solutions for equation 
(1.2) without the convective term ud x u is known for any p > 0, see e.g. [10, 11, 24, 27]. 
However, based on the reaction-diffusion experience mentioned above, one may expect 
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that the convective term prevents blow-up here as well. We will show below that this is 
indeed the case. Namely, the following theorem can be considered as the main result of 
the paper. 

Theorem 1.1. Let < p < |. Then, for every initial data uq G L 2 {VL) with zero mean, 
problem (1.2) possesses a unique solution u(t) which exists for all t > and remains 
bounded when t — > oo. 

Let p > 2. Then, there are initial data uo G L 2 (Q) with zero mean such that the 
corresponding solution blows up in finite time. 

Note that, in contrast to the situation with the semilinear heat equations, the result 
of Theorem 1.1 is not complete in the sense that we know nothing about the behavior 
of solutions for | < p < 2. Indeed, the sharp result for the semilinear heat equations 
stated above is strongly based on the maximum principle which we do not have for the 
CH equations, so we have to use alternative less powerful methods. We hope to return to 
this problem somewhere else. 

The paper is organized as follows. 

In Section 2, we introduce the notations and main technical tools which are necessary 
for our proof of Theorem 1.1. 

Section 3 is devoted to the proof of the blow up for the case p > 2. Actually, in the 
non-convective case, sufficient conditions of blow up of solutions for the CH equation can 
be established by concavity method of Levine (see [21]) rewriting the equation in the form 
of a nonlinear differential operator equations of the form 

Pu t + Au = F(u) (1.3) 

in a Hilbert space (see [16]). Here P, A are positive self-adjoint operators and F(-) is a 
nonlinear gradient operator . But the convective CH equation (1.2) can not be written in 
the form (1.3), so the method does not work directly and its adaptation to our situation 
requires some more delicate arguments and works only under the extra assumption p > 2. 

In Section 4, we prove the global existence and dissipativity of solutions of (1.2) in 
the case < p < g. Our proof uses the so-called Goodmann trick (see [15]) which is 
now-a-days the standard (for the theory of Kuramoto-Sivashinki equation) method to 
"extract" the dissipation from the convective term. However, in order to compensate the 
concave term w|-u| p , we need to chose the auxiliary function in this method depending on 
the size of the initial data and then use the so-called Gronwall lemma with parameter (see 
[13] and also [33], where this lemma was implicitly used to establish the global existence 
of solutions for the Navier-Stokes equations in a strip). Note that our approach also 
essentially employs the improved construction of the auxiliary function given in [3]. 

Finally, in Section 5, we consider some related problems which can be straightforwardly 
solved using the methods developed in the paper. In particular, we study here the problem 
of obtaining uniform in 5 — > upper bounds for the attractor of the so-called Kolmogorov- 
Sivashinsky-Spiegel equation 

d t u + d*u + d 2 x (2u-5u 3 )+ud x u = 0. (1.4) 

Furthermore, the so-called sixth order convective Cahn-Hilliard equation 

d t u - d*(d 2 x u + u - u 3 ) +ud x u = 0. (1.5) 
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are considered there. 

2. Notations and preliminaries 

In this section, we introduce the notations which will be used throughout the paper 
and introduce some technical tools important for what follows. 

As usual, we denote by H m = H m (—L,L) the Sobolev space of distributions whose 
derivatives up to order m belong to L 2 (—L,L). We write H = H(—L,L) instead of 
H° = L 2 (—L, L) and (•, •) stands for the usual scalar product in the Hilbert space H . 

The closure of C °°(-L,L) in H m (—L, L) will be denoted by H™ = H™(-L,L) and 
H™ r = H™ r (—L, L) stands for the subspace of H m which consists of 2L periodic functions. 
This definition works only for m G N, for negative or/and fractional ms, the correspond- 
ing Sobolev spaces are defined in a standard way using the duality and interpolation 
arguments respectively, see e.g. [31]. 

For every u G H™ r , m G N, we introduce the mean value operator 



1 f L 
( u ) : =2Lj L u ( x ) dx 



and denote by H™ r = H™ r (—L, L) the subspace of H™ r which consists of functions with 
zero mean: 

k"; er := {u G H£., (u) = 0}. 

We also introduce the inverse of the Laplace operator P := (— d 2 )' 1 defined on the 
functions with zero mean. It is well known that this operator gives the isomorphism 
between the spaces H™ r and H™+ 2 for every m G N: 

p . prm . rVm+2 p rrm _ Tjm+2 

' per per i per per 

and the following relations hold for all u G H 1 : 



\\u\\h < d \\d x u\\ H , \\P»(d x u)\\ H = \\u\\h, d := -. (2.1) 

7T 

Applying the operator P to the both sides of equation (1.2), we rewrite it in the equivalent, 
but more convenient (for our purposes) form 

Pd t u - d 2 x u - P(ud x u) = u\u\ p - (u\u\ p ). (2.2) 

We say that a function u(t,x) is a weak solution of (1.2) (or equivalently of (2.2)) on the 
time interval t G [0, T] if 

u G C([0,T],H per ) n L 2 ([0,T},H 2 er ), u G L p+1 ([0, T], L p+1 (—L, L)) (2.3) 

and equation (2.2) is satisfied in the sense of distributions. 

The following theorem gives the local well-posedness of this problem. 

Theorem 2.1. Let < p < 4. Then, for every u G H per , the problem (2.2) possesses a 
unique weak solution defined on the time interval t G [0, T], where T = T(uq) > depends 
on the H-norm of the initial data u$ 

The proof of this theorem is standard and, by this reason, is omitted. 
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Remark 2.2. As follows from the embedding theorem, 

ueC([0,T],H per )nL 2 ([0,T],H 2 per ) 

implies that u G L 10 ([0, T], L 10 ). Therefore, the assumption p < 4 guarantees that the 
nonlinear term u\u\ p belongs to L 2 ([0, T], H) and therefore it is subordinated to the linear 
terms in the equation. 

Let us mention also that the usual parabolic smoothing property works for such weak 
local solutions, so the factual smoothness of the solution u(t) for t > is restricted only by 
the smoothness of the non-linearity u\u\ p at u — 0. In particular, since this nonlinearity 
is of at least C 1+a for some a > (depending on p), one can show that we have at least 
u(t) G C 3 (— L,L) for t G (0,T] and any weak solution constructed in Theorem 2.1. Since 
this regularity is more than enough to justify all estimates used in the paper, we will not 
return to the questions of local well-posedness in what follows and will only concentrate 
ourselves on derivation of the a priori estimates which guarantee the global well-posedness 
or finite-time blow up of solutions. 

In the case p > 4, one can obtain the similar local well-posedness result just using more 
regular solutions, say, 

«GC([0,T],ij; er )nL 2 ([0,T],ij p 3 er ) 
and starting from more regular initial data u G H^ er . 

We conclude this section by stating two crucial lemmas, one of them will allow us to 
find sufficient conditions for blow up of solutions of the CH equation (1.2) and the other 
one will give the part of Theorem 1.1 related with the global solvability. 



Lemma 2.3. ([2.1]) Let ty(t) be twice continuously differentiable function that satisfies 
the inequality 

*"(*)*(*) - (1 + a) [V(t)} 2 > 0, t > 0, (2.4) 

and 

tf(0) > 0,tf'(0) > 0, (2.5) 
where a > is a given number. Then there exists 

such that 

ty(t) ->■ oo as t ->■ t\. 

Lemma 2.4. ^[13, 25] j Suppose that a > (3 > 1 and 7 > are given numbers that satisfy 
the inequality, 

^ < -!-, (2.6) 

a-l 7 + 1' V ; 

Suppose also that ^ is a non-negative absolutely continuous function on [0, 00) which 
satisfies, for some numbers K > 0, e Q > 0, M > and for every e G (0, £ ) the differential 
inequality 

+ e*(t) < Ke a \V(t)f + Me^ (2.7) 
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Then there exist a monotone function Q : IR + — > M + and a positive number k > such 
that 

< Q(^(0))e- K ' + Q(M). (2.8) 

Moreover, the dissipative estimate (2.8) remains true if the function \& > is only con- 
tinuous and satisfies the integrated version of inequality (2.7) 

< ^(0)e- £ * + Ke a [ e- £( *- s) [^(s)] /3 rfs + Me^- 1 (2.9) 

/or a// 1 > 0. 

3. Blow up of solutions to convective CH equations 
In this section we consider the following problem in Q := [-L, L\. 

d t u + d 2 x {p 2 x u + u 3 ) + ud x u = 0, 



u\ t=0 = u . 



(A) 



endowed by periodic boundary conditions. According to Theorem 2.1, this problem has 
a unique local weak solution for all uq G H per and it is equivalent to the following one: 

[Pd t u-d 2 x u + P(ud x u) =u 3 - (u 3 ), 

see Section 3. 

Our aim is to show that for some class of initial functions the corresponding solutions 
blow up in a finite time. For simplicity, we restrict ourselves to consider only the case 
p = 3 in equation (1.2), although as it is not difficult to see, the similar arguments work 
for all p > 3. 

The main result of the section is the following theorem. 

Theorem 3.1. Suppose that u is a solution of the problem (Ai) corresponding to the 
initial data u Q e H per which is not equal zero identically and 

Eo ■= ~\\P 1/2 u \\ 2 - l\\d x u \\ 2 + 1) > 0, (3.1) 

where A = \ + "&d\ and d := max{J , 1}; see (2.1). 
Then there exists t 1 < oo such that 

\\P l,2 u(t)\\ ^ oo, as 
Proof. We make the change u = e xt v. Then the function v solves the problem 
Pd t v + XPv - d 2 x v + e xt P(vd x v) = e 2Xt (v 3 - (v 3 )), 



\t=o= u °- 



(3.2) 



Let us consider the function 



*(*):= f \\P l l 2 v{r)\\ 2 dr + C Q) (3.3) 
Jo 
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where C is some positive parameter to be chosen below. Clearly 

tf'(i) = || J R 1 / 2 -^(t) || 2 = 2 [\pd T v,v)dT+\\P 1/2 u \\ 2 . (3.4) 

Jo 

Employing the equation (3.2) we also obtain that 

V"(t)=2(PdtV,v) = -2A||P 1/2 t;|| 2 -2||^|| 2 -2e A *(P(^),t;) + 2e 2At (t; 4 ,l). (3.5) 
By using the Cauchy inequality with e and the Schwarz inequality we obtain 

2e xt \(P(vd x v),v)\ = 2e xt \(P 1/2 (vd x v),P 1 / 2 v)\ < 2e xt \\pV 2 (yd x v) \\ \\P 1/2 v\\ < 

< 2d e xt \\v 2 \\ WP^vW < e ie 2Xt (v\ 1) + ^P 1 ^] 2 . 

Thus (3.5) implies 

> - (2\ + \\P 1/2 v\\ 2 - 2\\d x v\\ 2 + (2 - £i)e 2A V, 1). (3.6) 

Multiplying the equation (3.2) by d t v and integrating over (-L, L) we obtain the second 
main energy equality : 

jE{t) = \\P^d tV f + ±e 2Xt (v\l) + e xt (P(vd x v),d t v), (3.7) 

where 

E(t) :=-^\\P l ' 2 v\\ 2 -\\\d x v\\ 2 + \e 2X \v\l). 
Let us estimate the last term on the right hand side of (3.7) : 

e xt \(P(vd x v),d t v)\ <e At d |k 2 ||||P 1/2 ^|| < e 2 \\P 1 l 2 d t vf + ^e 2Xt {v\l). (3.8) 

Ae 2 

If 

A > <L (3.9) 

then it follows from (3.7) and (3.8) that 



f t E(t)>(l-e 2 )\\P^d t v\\ 2 . 



Hence 



E(t) > (l-e 2 ) f \\P l l 2 d T v{r)\\ 2 dr + E . (3.10) 
Jo 

It follows from (3.6) that 



*"(*) > 4(2 - e 1 )E(t) + 



2(i- ei )A-^ 



||P 1/2 H| 2 + (2-2 £l )||^| 



Let us take S\ — \ (it is important to have 4(2 — S\) > 4). Then we have the following 
estimate for fy"(t): 

*"(*) > 6P(t) + (A - 2d 2 )#'(t) + ll^f. (3.11) 
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Due to (3.10) and (3.11) with e 2 — | (note that for this e 2 the inequality (3.9) is 
satisfied) we have: 

> 5 r||P 1/2 ^(r)|| 2 rfr+(A-2^)^'(t). (3.12) 
Jo 

Here we have used the condition E(0) = E > 0. By using (3.4)and (3.12) we obtain the 
following inequality 

- ^ (^'(t)) 2 > 5 ^jf ||P 1/2 9 T t;|| 2 rfr^ ^jf \\P l ' 2 v\\ 2 dT + + 
- 5 Qf (P 1/2 8rU, P 1/2 u)<*r + ^||P 1/2 u || 2 ) + (A - 2dg)*tt' = 
jf ||P 1/2 9 T t;|| 2 rfr + (f) - (^j\p 1/2 d T v, P l ' 2 v)dT + 



pl/2,. 112^ 2 ' 
r Uq\\ 



2 

2\ 



+ 



(\-2dl)y(t)y'(t)-C V(t). (3.13) 
Let us choose in (3.13) 

C7o = ^||P 1/2 «oH 2 . 

Then the expression in the square brackets is nonnegative due to Cauchy-Schwatz in- 
equality. Therefore we have 

- \ (*'(0) 2 > ( A - 2dg)tf (*)*'(*) - C tt(t). (3.14) 
According to (3.11) the function ^f'(t) is a non-decreasing function. Thus 

tt'(t) = iip 1 /^^)!! 2 > ||p 1/2 m || 2 . 

Hence we obtain from (3.14) 

tt"(t)*(t) - ^ W(t)] 2 > (A - 2d 2 )M>(t)||P 1/ V)|| 2 - C *(f) > [2(A - 2d 2 ) - l]C„tt(f). 
Finally noting that A = | + 3d 2 , we obtain: 

*"(*)*(*) - \ Wit)] 2 > o. 

Hence due to the Lemma 2.3 the statement of the Theorem 3.1 holds true. □ 

Remark 3.2. It is clear that Theorem 3.1 is true for solutions of the equation (A) under 
the homogeneous Dirichlet's boundary conditions. 

Remark 3.3. It is easy to see that the result of the Theorem 3.1 remains true also for 
the multi-dimensional convective CH equation of the form 

d t u + A(Aw + u 3 ) + ub.Vu = 0, 

where b G L°°(Q,M.) is a given vector field and u satisfies the appropriate boundary 
conditions. 
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4. The convective CH equations: global existence 

In this section, we continue our study of the convective CH equation: 

d t u + ud x u + d 2 x {d 2 x u + u\u\ p ) = (4.1) 

on the interval Q = [—1,1] (for simplicity, we take L — 1 here) endowed with the periodic 
conditions. We have seen in the previous section that this equation possesses the blowing 
up in finite time solutions if p > 3. The aim of the present section is to show that the 
presence of the convective term prevents the blow up if the exponent p is not large. 
Namely, the following theorem is the main result of the section. 

Theorem 4.1. Let the exponent < p < |. Then, for every u G L 2 ([— 1, 1]) with zero 
mean, problem (4.1) possesses a unique solution defined for all t > and the following 
estimate holds: 

\\u(t)\\ L 2 < Q(\\u \\ L 2) e - at + C„ (4.2) 

where the positive constants C* and a and the monotone increasing function Q are inde- 
pendent oft and u Q . 

Proof. Since the local existence and uniqueness theorem for the equation (4.1) equation is 
standard and immediate, we only need to verify the dissipative estimate (4.2). As usual, 
we start with the case of odd periodic solutions and use the following Lemma which is in 
fact proved in [3] 

Lemma 4.2. For every sufficiently large N there exists a 2-periodic function <fi with zero 
mean such that 

WWh* < CN Z ' 2 , \\(t>\\ L oo <CN (4.3) 
with constant C independent of N, such that, for every u G H 2 (Q) with u(0) = 0, 

\\u xx \\ 2 L2 - ((j) x ,\u\ 2 ) > N\\u\\ 2 L2 . (4.4) 

The function with the desired properties is constructed up to scaling in [3]. Indeed, 
it is proved there that, for any sufficiently large L, there exists a 2L-periodic function 
V> G H 2 (-L,L) such that 

\\^\\ H2 < CL 3 ^ 2 , < CL 

and, for any u G H 2 (—L, L) with u(0) = 0, the following inequality is satisfied: 

1 f L 

i\d 2 u(x)} 2 — u 2 (x)d x if) , dx > - / u 2 (x)dx. 
-l 2 J_ L 

Scaling x = Ly, = L 3 ip and iV = CL 4 , we end up with (4.4). 

We now return to the key a priori estimate (4.2) for the odd periodic solutions of (4.1). 
To this end, for any large N, we multiply equation (4.1) by v :— u — </> where = 0at is 
constructed in Lemma 4.2. Then, after some transformations, we get 

dtWvfv + 2\\d 2 x u\\ 2 L2 - (d x <f>, \u\ 2 ) = -2(u\u\ p , 8 2 x u - 8 2 X <P) + 2(d 2 x u, d 2 x <f>). (4.5) 

Using the Cauchy-Schwartz inequality together with (4.3) and (4.4), we get 

dt\Hh + 1/2||^||| 2 + iV|M|| 2 < \\u\\ 2 <*£l + CN\ (4.6) 
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Note also that, due to Lemma 4.2, we see that, for every q e [1, oo], 

\\u\\ Lq -CN < \\v\\ Lq < \\u\\ Lq + CN, (4.7) 

where the constant C is independent of N. 
Applying the interpolation inequality 

to estimate the right-hand side of (4.6), we end up with 



dt\\v\\ 2 L 2 + iV||u|||a < C\\u\\ 2 L t P + CN 3 . (4.8) 

In order to derive the desired dissipative estimate for u from (4.8), we'll use the Lemma 
2.4. To this end, we scale time t = N 2 t and introduce e — 1/N. Then, (4.8) reads 



3p+4 



>3p+4 



d r\\v\\ 2 L 2 + e||w||i 2 < Ce 2 \\u\\ L t p + Ce" 1 . 

Integrating this inequality in time, introducing *&(t) := ||w(t)||^2 and using (4.7), we end 
up with 

tf(t) < ^(0)e- £ * + Cs 2 / e- £ ^- s \m{s)} W rf s + Ce~ 2 , (4.9) 
where C is independent of e — » 0. Equation (4.9) ha the form of (2.9) with a = 2, 

3p+4 
4— p 



(5 = and 7 = 1 and the condition (2.6) reads 



4p 1 

< 



4-p 2 

This condition is i satisfied if and only if p < |. Thus, due to Lemma 2.4, the odd solutions 
of (4.1) cannot blow up in finite time and satisfy the dissipative estimate (4.2) if p < |, 
so in the particular case of odd initial data, Theorem 4.1 is proved. 

We are now ready to consider the general case of periodic u G H per with using the so- 
called Goodman trick, see [15]. Namely, we consider a circle of shifted functions <f> s (x) : = 
<f)(x + s), s 6l and introduce 



R(t) := min \\u(t) — (j). 



|2 

5 1 • 



Then, due to (4.3), we have the analogue of (4.7): 

IK*)llL -CN < R(t) < \\u(t)\\ 2 L 2 + CN (4.10) 
and for the minimizer s (t), we have 

(u(t),d x cf> s{t) ) = 0, (4.11) 
and, at least formally (see [15] for the justification), 

\jt R ® = \ jt [][U ^ > ~ ^<*>H* a = ( d Mt),u(t) - <p s{t) )- 

- s'(t)(d x <i> 8 (t),u(t) - (j) s (t)) = (d t u(t),u(t) - s(t) ). 
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Thus, multiplying equation (4.1) by 2(u(t) — </> s (t)) an d arguing as before, we get 

jR{t) + 3/2||^||| 2 + (9^ s(t) , \u\ 2 ) < \\u\$$, + CN 3 (4.12) 

However, in contrast to the odd case, we cannot apply directly Lemma 4.2 since the 
condition u(t,s(t)) = is not necessarily satisfied. So, we need to introduce a time 
dependent "constant" c(t) := u(t,s(t)) and a function w(t) = u(t) — c(t) for which the 
conditions of the Lemma 4.2 are satisfied, and we have 

\\d 2 x w\\ 2 L i + (d x <f> s{t) ,\w\ 2 )>N\\wf. 

Note that, due to the zero mean condition on u, 

\\w\\ 2 = \\u\\ 2 L 2 + Mc(t) 2 > \\u\\ 2 

and, due to the orthogonality condition (4.11), 

{d x <t> s (t), M 2 ) = (dx(f> s (t), H 2 ) - 2c(t)(d x <t> a (t),u(t)) + c(t) 2 (d x (j) s{t) , 1) = (d x <p s{t) , \u\ 2 ). 
Therefore, 

\\dlu\\ 2 L , + (9 x s(t) , M 2 ) = \\d 2 x w\\ 2 L2 + (d x cf> s{t) , \w\ 2 ) > N\\u\\ 2 
and (4.12) implies that 

A 2 3p+4 

-R(t) + N\\u\\ 2 L2 < C\\u\\ L r + CN 3 . (4.13) 

Finally, integrating (4.13) in time, using (4.10) and arguing exactly as in the case of odd 
initial data, we derive the desired dissipative estimate (4.2) for general Uq G H per and 
finish the proof of the theorem. □ 

5. Related problems 

In this section, we apply the above considered methods to some equations which 
are, in a sense, close to the Kuramoto-Sivashinski and Cahn-Hilliard equations, such 
as Kolmogorov-Sivashinski-Spiegel equations and for the sixth order convective CH equa- 
tions. 

5.1. Convective CH Equation vs Kolmogorov-Sivashinsky-Spiegel equation. 

We consider now the problem 

'd t u + d*u + d 2 x (2u-5u 3 )+ud x u = 0, 
«lt=o = u o 

in the domain Q := (-L, L) endowed by the periodic boundary conditions. We are 
going to show that the estimate obtained in [12] for the size of the absorbing ball can be 
improved at least for small values of 5 where one expects that the Kuramato-Sivashinsky 
dynamics will dominate. Indeed, to the best of our knowledge, all previous methods of 
obtaining the dissipative estimates for this equation utilize only the dissipativity which 
comes from the cubic term ignoring the extra dissipation provided by the convective term. 
As a result, the obtained estimates were divergent as 5 — > 0, see [12, 6, 7, 8]. In particular, 
the radius of the absorbing ball in H per constructed in [6] behaves like L 3 5~ 1 / 2 as 8 — > 0. 
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Using below the technique related to the Kuramoto-Sivashinski equation (analogous to 
what is used in Section 4, we show that the radius of the absorbing ball remains bounded 
as 5 — > 0. Since we are not interested in the dependence of this radius on L, we set L = 1 
for simplicity. Then, the following theorem holds. 

Theorem 5.1. Let L — 1. Then, for every uq G H per , problem (4.1) possess a unique 
solution and the following estimate holds: 

\Wt)\\h < C\\u \\ 2 e- at + C* (5.2) 

where the positive constants C , a and C* are independent of 5 — > 0. 

Proof. The existence and uniqueness for that equation is well-known, see e.g., [8], so we 
only derive the uniform estimate (5.2). Analogously to Section 4, we start with the case 
of odd initial data u and multiply (5.1) by v = u — if where ip = tp^ is the same as in 
Lemma 4.2 and iV will be fixed below. Then, after the obvious transformations, we have 

dt\\v\\ 2 L 2 + 2\\d 2 x u\\ 2 L , - (d x <f>, \u\ 2 ) + Q5(u 2 , \d x u\ 2 ) = 

= Q5(u 2 d x u, d x <f>) + M\d x u\\ 2 L2 - 4(d x u, d xV ) + 2(d 2 x u, d 2 J). (5.3) 

We estimate the first term in the right-hand side of (5.3) via the Cauchy-Schwartz in- 
equality: 

65(u 2 d x u, d x <j>) = Q5(ud x u,ud x <p) < 65(u 2 , \d x u\ 2 ) + 65(|^0| 2 , u 2 ) < 

< 65(u 2 , \d x u\ 2 ) + CN 3 5\\u\\ 2 L2 , (5.4) 

where we have used that || d x (f>\\ < C||0||^2 < CN 3 ^ 2 . The rest terms in the right- 
hand side can be estimated in a standard way using the interpolation inequality < 

ClIHIfra IMIl^ an( i the Cauchy-Schwartz inequality. Therefore, using again (4.3), we see 
that for sufficiently large N, 

dt\\v\\ 2 L 2 + (f3N - a{\ + 5N 3 ))\\u\\ 2 L2 < CN 3 . (5.5) 

where the positive constants j3, a and C are independent of 5 and N. Finally, fixing 
such that (5N > 3a, we see that, for 5 being small enough that 5N 3 < 1, the following 
inequality holds: 

d t \\v\\ 2 L 2+a\\u\\ 2 L 2 <CN 3 

and the Gronwall inequality applied to it gives the desired uniform estimate (5.2). Thus, 
the theorem is proved in the particular case of odd initial data. The general case can be 
reduced to that particular one using the Goodman trick again, exactly as in Section 4. 
So, the theorem is proved. □ 

Remark 5.2. It looks natural to consider the mixture of problems (4.1) and (5.1), namely, 

d t u + uu x + (u xx + u\u\ p - 5u 3 ) xx = 

with p < | and study the limit <5 — ?► 0. However, the above method does not work at 
least directly in this case (at least without stronger assumptions on p). Indeed, as we see 
from (5.5), if we allow N to be dependent on the norm of the initial data (as in Section 
4), we also need to decrease 5 in the dependence on the initial data. Thus, the problem 
of obtaining the uniform with respect to 8 — > estimates for that case remains open. 
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5.2. Sixth order convective Cahn - Hilliard equation. We consider here the follow- 
ing sixth order convective Cahn-Hilliard equation 



on a bounded interval fl = (-L, L). This equation was derived in [29] as a model of process 
of growing crystalline surface with small slopes that undergoes faceting where u = d x h 
is the slope of a surface h(x,t). In [18] the authors proved existence and uniqueness of 
the global solution to initial boundary value problem for (5.6) under periodic boundary 
conditions. 

We are going to show that the semigroup generated by the initial boundary value problem 
for (5.6) under the periodic boundary conditions is dissipative (the case of, say, Dirichlet 
boundary conditions can be treated similarly), namely, the following proposition holds. 

Proposition 5.3. For any u e H~^ r , problem (5.6) possesses a unique solution u e 
C(R+, Hper) n LfoJ^+i Hper) an d ^ e following dissipative estimate holds: 



for the positive constants a and C* and monotone function Q which are independent of t 
and uo- 

Proof. Applying to both sides of (5.6) the operator P = (—dl)' 1 we get an equivalent 




(5.6) 




(5.7) 



problem 




(5.8) 



Multiplying the equation (5.8) by u and integrating over (0, L) we obtain 




Employing the inequality (2.1) and the interpolation inequality 





P^u\\ 2 + (1 - £i)||c^|| 2 + 3(m 2 , (d x u) 2 ) < d ||u 2 ||||P3 U || + CJP^II 2 < 




. 2 || M 2 || 2 + (d + C7 2 )||P*u|| 2 < 25 2 || M 2 || 2 + [MC \ + ° 2)] 



L. (5.9) 




P h u f + (i _ £l )||^ u || 2 + (3 - 8e 2 d 2 ) (u 2 , (d x u) 2 ) < K . 



(5.10) 



CONVECTIVE CAHN-HILLIARD EQUATION 



13 



By choosing e± — \ and e 2 = and applying the Gronwal inequality, we deduce (5.7) 
and finish the proof of the proposition. □ 

Remark 5.4. Using the standard parabolic regularity, one may show that the solution 
u(i) constructed in Proposition 5.3 becomes C°° (and even Gevrey) regular for all t > 0. 

Remark 5.5. By using the same arguments as in the proof of the Theorem 3.1 we can 
show that a wide class of solutions to the initial boundary value problem for the sixth 
order convective CH equations with concave potential 

d t u-d*(d 2 x u + u + u 3 )+ud x u = 0. (5.11) 

blow up in a finite time. 

We would like also note that the blow up theorem for sixth order unstable CH equations 
(5.11) without the convective term can be established by using the concavity method of 
Levine since in this situation the equation can be written in the form (1.3). 

Finally, arguing as in the proof of Theorem 4.1, one can show the global existence and 
dissipativity of solutions of the following problem: 

d t u - dl(d 2 x u + u + u\u\ p ) + ud x u = (5.12) 

if p < po for some exponent po > 0, so the presence of the convective term prevents blow 
up in that situation as well. However, in order to compute the exponent po, we need the 
analog of the sharp Lemma 4.2 for the six order operator which we do not present here. 
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